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ABSTRACT 



A Bayesian analysis was used to determine the host star parameters of 2740 
candidate and confirmed Kepler transiting planets and estimate the probability 
Phz that a planet lies within its star's habitable zone, i.e., that the incident stellar 
irradiation lies between the threshholds for the runaway wet greenhouse and polar 
CO 2- condensing "snowball" climate states. This analysis is based on photometry 
from the Kepler Input Catalog, an ensemble of Dartmouth Stellar Evolution 
Program models, and priors for stellar mass, age, metallicity, distance, as well 
as a constraint from the duration of the planet transit, and yields probability 
density functions for stellar parameters with which confidence intervals of planet 
radius and incident stellar flux as well as pnz are calculated. Fifty- two planets 
have Phz > 0.5 and a most probable stellar irradiation within the habitable zone 
limits. Five of these have estimated radii less than twice the Earth, but the 
majority are super-Earth to Neptune size, and 8 are gas giants which could host 
habitable exomoons. Another is possibly not a planet, but a brown dwarf or 
ultracool star. The planet candidate most resembling Earth in terms of size and 
irradiation is KOI 3010.01, which orbits a late K or early M dwarf. Parallaxes 
from the Gaia mission will reduce errors in planet radius and irradiance by at 
least a factor of 5, permitting definitive assignments of transiting planets to 
habitable zones of Kepler stars and, eventually, studies of any correlations of 
planet properties with position in those zones. 
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Subject headings: planetary systems 



— techniques: transit surveys 



— astrobiology 
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Introduction 



The Kepler mission was launched in March 2009 with a mission to find Earth-siz e 



plane ts in the circumstellar "habitable zone" (HZ) of solar-type stars f lBorucki et al 



20101 ) ■ Broadly speaking, the HZ is considered the range of orbital semimajor axes over 
which the surface temperature on an Earth-like planet would permit liquid water. A 
more narrow definition is that it is the range of stellar incident fiux between the runaway 
"wet" greenhouse limit - beyond which a water-vapor saturated N2-CO2 atmosphere 
cannot radiate, and the CO2 condensation limit below which atmospheric CO2 - the 

larth-like atmsophere - condenses at the poles 



primary green 


rouse 


(Kastine et al. 


1993 



Ishiwatari et al. 



20071 ). This de finition is subject to assumptions about 



20 planeta ry albedo, orbita l eccentricity and obliquity 



oceans (Abe et al. 



WiUiams fc PoUardl fl2003h. extent o f 



201 ll ). and composition of the atmosphere (jPierrehumbert et al. 



20111 ). 



Borucki et al. 



(120111 ). assuming an Earth-like Bond albedo of 0.3, published a catalog 



23 of 54 (out of 1235) candidate planets or Kepler Obiects of Interest (KO I s) wit h equilibrium 



24 emitting temperatures between 273 and 373 K. 



Kaltenegger &: Sasselovl (1201 ll ). noting the 



25 importance of albedo, specifically cloud cover, to equilibrium temperature, co mputed inn er 



26 and o uter HZ boundaries based on the stellar irradiation criteria derived by 



Selsis et al. 



27 (I2OO7I ) for high H2O and high CO2 atmospheres, respectively. They identified 76 possible 



28 habita ble planets, de p endin g on the assumed fractional cloud cover. They found that many 



29 of the 



Borucki et al. 



( I2OIII ) candidates were too hot for this habitability criterion and 

30 pointed out that errors in stellar parameters contribute most to the uncertainty of whether 

31 a planet orbits within the HZ. 

32 Subseque ntly, a larger sample of 2300 KOIs, including some that have been confirmed, 



33 was released ( Batalha et al. 



20121 ). Stellar parameters for KOI host s, i.e., mass M^, and 



34 radius i?*, were determined by fitting Yale-Yonsei model isochrones ( iDemarque et al 



200J) 
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35 to values of effective temperature (T*), surface gravity (log (7), and metallicity ([Fe/H]). The 

36 latter were derived from the photometry of the Kepler Input C atalo g (KIC) and a r nodel of 



37 stellar populations and Galactic structure (IBrown et al 



201l[ ). The 



Batalha et al. 



(120121 ) 



38 estimates of mass and radii assumed gaussian-distributed errors and employed standard 

39 deviations derived from a comparison b etwee n KIC-derived parameters and spectroscopic 



40 values. They r evised the iBrown et al. 



Batalha et al. 



(120111 ) estimates of log (7 and -R* for many stars. 



(|2012[ ). assuming an albedo of 0.3 and full redistribution of heat over a 



42 planet's surface, identified 46 candidates with 185K < Te„ < 303K 



However, the 



Brown et al. 



(120111 ) stellar parameters themselves are uncertain and, 

44 in some aspects, problematic. The vast majority of Kepler stars do not have measured 

45 parallaxes, at least for now. KIC-derived effective temperatures are about 200 K hotter than 

46 estimates based on the i nfrared flux method and K IC photometry must be corrected to place 



47 it in the Sloan system ( iPinsonneault et al 



2OI2I ). Moreover, errors in stellar parameters 



4B can be strongly non-gaussian, particularly for solar-type stars for which ph otometry is 
49 unab l e to distinguish 



2011 



Gaidos fc Mann 



jet we en main sequence and evolved (subgiant) stars (iBrown et al. 



20131 ). In such cases, standard deviations have limited utility. 



51 A revised and expanded catalog of 2740 confirmed and candidate planets around 2036 

52 Kepler stars was released on 7 January 2012. Here, I identify 52 candidate HZ planets 

53 by comparing corrected KI C qrizJHK photorn etry w ith isochrones from the Dartmouth 



54 Stellar Evolution database ( iDotter et al 



20081 ). Like 



Brown et al 



55 statistics, but with new priors that describe dist ributions with mass 



56 and distance using recent models of the Galaxy (jVanhoUebeke et al. 



(l201l[ ). I use Bayesian 
IMF ), metallicity, age. 



20091 ). In addition, the 



57 duration and probability of the planet transit itself are used to constrain the radius of the 



58 host star ( Plavchan et al. 



2OI2I ). This statistical approach is used to calculate a probability 



59 that any given planet is in the HZ. A similar approach has been used on Hipparcos stars 
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( iBailer- Jones 



201l[ ). Alternative a pproaches to the use o f bro ad-band photome t ry to derive 



61 stellar parameters are described in 



Ammons et al. 



(120061) and 



Belikov fc Roserl fcoOSf ) 



Methods 



2.1. Algorithm 



1 compare photometry for each star with sets of synthetic SDSS qrizJH K,^ photometr y 



65 from the isochrones of the Dartmouth Stellar Evolution Database (DSED) ( iDotter et al. 



20081 ). Specifically I compare up to six colors constructed with respect to the r magnitude. 



According to Bayes theorem, the probability Pi that the ith model (hypothesis) is supported 
by the photometry is equal to the probability that the colors cj can be produced by the 
model, multipled by a prior function pi. Assuming gaussian-distributed errors in photometry 
that probability is, 



2a] 



(1) 



where the summation is over up to six colors, cj are the synthetic colors, aj are photometric 
errors, kj is the interstellar reddening coefficient for the color, and is the amount of 

reddening that is assigned to a particular model and star (see below). The normalization 
in Eqn. [T]is unimportant as it independent of the models and 1 seek the model which has 
the largest value of P. The prior is the product of individual priors for the mass, age, 
distance, and metallicity of the model, the intervening extinction, and, since at least one 
planet has been detected around each of the se stars, a constraint on stellar radius placed 



by the duration and probability of a transit (jPlavchan et al. 



2012h . 



Photometry and other data for the host stars of the KOIs were extracted from the 
KIC catalog available at the MAST database. KI C qriz magnitudes were co rrected to the 



81 Sloan system using the corrections determined by 



Pinsonneault et al. 



mm . Standard 
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errors for each bandpass were estimated using the expression a = aolO^'^ mo)/2.5^ where 



(To = 0.02, 0.02, 0.015 , 0.015, 0.02,0 



84 rizjHK, respectively (IBrown et al 



025,0.02, and mn = 15 , 15, 15.3, 15.3, 13, 11.75, 10.8, for 



2011 



Cutri et al 



20031 ). For griz the magnitude m is the 



85 Kepler magnitude Kp and for JHKs it is the respective magnitudes. Errors in color were 

86 calculated by assuming that errors in individual bandpasses are uncorrelated and adding 

87 the two corresponding such errors in quadrature. 

88 Minimization of P with respect to Eb-v leads to a formula for the best-fit reddening 

89 for each model: 



E 



B-V 



(2) 



I adopted extinction coefficie nts A of 3.758 fg), 2 565 f r), 1.874 (i) 1.37 7 {z), 0.272 (J) 



0.173 (if), and (i^), based on 



Girardi et al. 



(120051) and 



Chen et al 



fl2007h . 



92 I used the DSED interpolator tool to construct a grid of isochrones with [Fe/H] 

93 G [-1.5, +0.5], spaced by 0.1 dex, a/Fe G [-0.2, +0.4], spaced by 0.2, and ages G [1,12] 

94 Gyr. All models used a helium fraction Y = 0.245 + 1.5Z, where Z is the total heavy 

95 element abundance. I further restricted the selection to stars with initial masses between 

96 0.1 and 2 solar m asses, as late M and O, B, and early A-type stars are absent from the 



Kepler target list flBatalha et al. 



2010l ). This restriction reduced the total number of models 



98 considered to 657,347. 

99 Once the best-fit model with the maximum P was found, additional models (typically 

100 a few dozen) with neighboring (difference less than 1.5 times the grid spacing) values of 

101 mass, age, [Fe/H], [a/Fe] were identified. A set of 100 linear interpolations between the 

102 best fit model and each of these neighboring models was made and probablities calculated 

103 using Eqn. [TJ The interpolation yielding the highest value of P was recorded. 
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104 2.2. Prior functions 

105 Priors essentially weight each DSED model, i.e. each combination of initial mass, 

106 metallicity, and age. No prior is imposed on the allowed range of a/Fe between -0.2 and 

107 +0.4. The distance modulus for each star/combination is computed in the r-band, i.e. 
lOB yU^. = r — ArEs-v — Mr . As a p r ior fo r initial stellar masses I adopt the tripartate initial 



109 mass function (IMF) of iKroupal (120021 ). Priors for age, metallicity, and distance modulus 
no fi = m — M were constructed using the distributions of dwarf stars (log (7 > 4) with Kp < 16 
111 synthesized using TRILEGAL . The simulated population was rest ricted to dwarfs to reflect 



the criteria of the selection of Kepler targets (IBatalha et al 



values for key TRILEGAL parameters are the same as used in 



2010h. The fmostlv defa ult) 



Gaidos fc MannI fl2013h . 



114 The resulting distributions (Fig. [T]) have a median metallicity of -0.13, median age 

115 of 3.9 Gyr, and median /i = 11.2 (~1740 pc). The age distribution is complex because 

116 the population includes halo stars, which form 11-12 Gyr ago, and disk stars, which start 

117 forming 9 Gyr in these simulations. The star formation rate in the disk is modeled as 

118 two steps, with the second occurring at about the epoch of the Sun's formation. The 

119 paucity of stars younger than 1 Gyr is partly due to the fact that the Kepler fleld probes 

120 the stellar population that is > 100 pc above the Galactic plane. Of course, stellar ages, 

121 metallicities, and distances are interrelated, but here they are used separately, providing 

122 broad constraints on the possible ranges of stellar parameters. The distance distribution 

123 is particularly important in allowing the flnite scale height of the Galactic disk to prevent 

124 Malmquist bias from selecting arbitrarily distant and luminous stars. 

125 The Kepler fleld is ~10 degrees wide and close to the Galactic plane {b ~ 13 deg), 

126 so the stellar populations that are probed will vary signiflcantly across the fleld. Using 

127 TRILEGAL, I synthesized the stellar population over a square degree centered at each of 

128 the 85 Kepler CCD flelds. Only those synthetic populations for CCD fleld centers with b 
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129 within 0.5 deg of a given Kepler star (about 10% of the total) were used to calculate priors 

130 for age, metallicity, and distance. 

131 A prior for extinction Eb-v necessarily involves information about both the distribution 

132 of stellar distances, and the distribuiton of dust along the line of sight. However, by 

133 assuming that the spatial distributions of stars and dust along the line of sight is the same. 



134 the prior becomes particularly simple: a uniform distribution bet ween and tota . 



Schlegel et al. 



(oo' 



fll998h 



135 extinction along the line of sight (see Appendix A). I adopted the 

136 Galactic reddening maps and interpolated the total extinction at the coordinates each star 

137 using the IDL tools provided by the Princeton website. Models with optimal Eb-v values 
13B outside this range are not disallowed, but reddening is limited to — Eb-v{.'^) and the 

139 models are penalized for the resulting disagrement between measured and model colors 

140 (Eqn. H]). 



2.3. Constraints from the planet transit 



The transit duration r a nd orbital period Pk of a transiting planet constrain stellar 



143 density ( iPlavchan et al 



20121 1 and can be used as an additional prior for stellar models. In 

144 the case of Kepler low-cadence data, the constraint is weakened by a lack of information 

145 about the orbit, specifically independent determination of the orbital eccentricity and the 



146 transit impact parameter. The transit duration D is: 



K 



1 + e cos < 



(3) 



147 where the stellar free-fall time is t = 2y R^, / (ttGM*), G is the gravitational constant, e is 

148 the orbital eccentricity, (p is the argument of periastron relative to the line of sight to the 

149 star, and b is the impact parameter. Given D and Pk and a value for r for each stellar 
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model, e can be written ClS db function of b and 



e(6,0,A) 



1-62) (l-62_A2si 



sm 



cos ( 



1 - 62 + A2 C0S2 



(4) 



where A = D / [j'^^^P^/^). The eccentricity was calculated over a uniform grid of 6 G [0, 1] 
and G [0,27r]. Each value of e was assigned a probability, i.e. values not G [0,1] are 
assigned zero and others are assigned probabilities from a prior distribution of e. A Rayleigh 
distribution, 



P(e) 



-£2/(2^2) 



(5) 



155 is assumed . Such a distribution has been used in a previous analysis of Kepler transit 



156 durations ( iMoorhead et al 



201 ll ) and is motivated by dynamical theory (IJuric fc Tremaine 



20081 ). Then the prior for the ith model from the duration of the transit is 

2tt 



Pi 



1 



db 



27r ./o JO 



d<f)p{e{b,<j);Ai)). 



(6) 



158 To account for finite errors in transit duration, the prior can be calculated using multiple 

159 Monte Carlo-generated values of D is repeated and then averaged. In the case of a 

160 multi-planet system j = 1...N, the product of the individual transit duration priors 



was used. A va 


ue of 


(Moorhead et al. 


2011) 



ue of a = 0.2/ ^yn/ 2 for the dispersion in eccentricities was used 



20111). Figure [2] plots p for 4 values of cr. The choice of prior 
for e, although motivated by Kepler observations, is not unique and indeed a more refined 
distribu tion would include d ependencies on planet mass, orbital period, and the age of the 



165 system (jwang fc ZhoJ 



201ll). 
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166 2.4. Probability that a planet orbits in the habitable zone 

167 Orbit-averaged irradiation is only weakly dependent on eccentricity for near-circular 

168 orbits. 1 assume near-circular orbits in which case the mean irradiation in terrestrial units is 

4/3 / \ 2/3 

Lq V 365.24 d) \M(r, 



I^!^(^JL_X {^\ . (7) 



169 A planet is defined to be in the HZ if lout < I < hni where the irradiance of th e inner edge 



170 of th e HZ for a 50% cloud-covered planet with efficient heat re-distribution is (ISelsis et al. 



20071) 



hn = [0.68 - 2.7619 X lO^^B - 3.8095 x 10^*^6^] ^ , (8) 

172 and the outer edge is: 

= [1.95 _ 1.3786 X 10-^6 - 1.4286 x lO'^G^] , (9) 

173 where 9 = T* — 5700. 

174 1 determine whether a planet is in the HZ for each set of model stellar parameters 

175 M*,L*, and with associated probability Pi. The probability pnz that the planet is in the 

176 HZ is then: 



Phz = ^ — (10) 

i 

177 I consider candidate habitable planets (or planets that may have habitable satellites) as 
17B having pnz > 0.5 and having a most probable value of / (with highest P) satisfying 

179 lout < I < lin- 
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Results 



3.1. Comparison with spectroscopic parameters 

Accurate estimates of stellar effective temperature T* and radius i?* are crucial to 
assessing whether a planet is in the HZ, as together these largely determine the luminosity 
of the host star and the irradiance experienced by the planet on a given orbit. The inferred 
radius of a transiting planet also scales linearly with the estimated radius of the host star. 
The radii of distant Kepler stars cannot be directly measured, but spectroscopic values of 



187 T^: and surfac e gravity log o, the 



188 with planets (IBruntt et al 



2012 



atter related to -R^,, a re available for some Kepler stars 



Buchhave et al 



2012 



; Mann et al., in prep.). Figs. |3] 



and m compare photometry-based values of and log g with reported spectroscopic values. 
Photometric values for solar- type st ars where all 6 colo r s are available average 210 K higher 
than spectroscopic values (Fig. [3]). 



Pinsonneault et al. 



( I2OI2I ) found that both the original 



KIC temperatures and spectroscopic estimates were ~215 K cooler than determinations 
using the infrared flux method (IRFM). Thus the new photometric estimates are in line 
with IRFM values. The offset between photometric and spectroscopic temperatures is less 
(~100 K) for M dwarfs; spectroscopic temperatures for these stars (Mann et al., in prep.) 
were det ermined by compa ring spectra to synthetic spectra from PHOENIX/BT-SETTL 



197 models flAUarc 



Boyajian et al. 



et al. 



f l2012h . 



20111 ) and tuning the comparison using the temperature estimates of 



After removing the 210 K offset, the standard deviation between photometric and 
200 spectroscopic val ues of for solar-ty pe stars is a = 205 K. This equals the performance 



201 of the analysis of iBailer- JonesI (l201lf ). but without the benefit of parallax. Photometric 

202 T^, for several stars is significantly (> 2a) lower than spectroscopic estimates. All but one 

203 of these are missing either i- or 2;-band photometry, or both. The importance of these 

204 bandpasses is not surprising as they are the only source of information in the wavelength 
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205 range 0.7/im < A < 1.1/im, just beyond the peak in emission from most of these stars, a 

206 spectral feature which most strongly constrains T*. 

207 There are also 7 stars which have significantly (> 2cr) hotter photometric estimates of 
20B T^,; only one of these has missing photometry. One of these has a spectroscopic temperature 

209 estimated using the method applied to M dwarfs (Mann et al., in prep.) but does not have 

210 the spectrum of an M dwarf, and thus its spectroscopic temperature is unreliable. The 

211 reason(s) for the discrepancy among the remaining six are unclear. One possibility is that 

212 the photometric source is a blend resolved by spectroscopy, or that the transit signal itself 

213 may be coming from a component of a blend which is dissimilar to the source of most of the 

214 light, and consequently the transit duration prior is skewing the stellar parameters. 

215 Photometry-based estimates of logg are more discrepant with spectroscopic values, 

216 although an overall correlation is apparent (Fig. H]). About half of the most discrepant cases 

217 lack photometry in at least one bandpass, although many stars with missing photometry 
21B are assigned surface gravities close to the spec troscopi c estimates. Althou gh photometric 



219 colors involving the SDSS u ( iLenz et al 



19981 ) and z (IVickers et al. 



20121 ) bands can be 



220 used to discriminate between hotter main sequence and evolved stars, photometry is a much 

221 blunter tool to separate solar-type stars by luminosity class. While my analysis may only 

222 marginally improve this situation, it does quantify the uncertainties. 



Among the KOI host star with reported 



spectroscpic pa r amete rs are those with 



Buchhave et al. 



( 120121 ) report spectroscopic 



224 candidate HZ planets discussed below ( §3.2p . 

225 parameters for 3 of the stars in Table [H including Kepler-22h. The photometric values 

226 of are up to 200 K warmer than the corresponding spectroscopic estimates (Fig. [3]). 



Muirhead et al. 



228 and Mann et a. 



Muirhead et al 



(120121 ) obtained infrared (i^-band) spectra for 6 of these HZ candidates 
(in p rep.) obtained visible-wavelength spectra for 10, including 4 of the 



( 120121 ) targets (Tabled]). The spectra show that all 12 are late K- or early 
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230 M-typ e dwarfs. Two (KOIs 812.03 and 854.01) were previously identified by 



Muirfiead et al. 



231 ( 120121 ) as orbiting in the star's habitable zone. In general, the photometric temperatures are 



232 ~100 K warmer than the spect roscopic values of Mann et al, which are in turn warmer than 



233 those of 



Muirhead et al. 



(120121 ) (see also Fig. [3]). The most discrepant case is KOI 2418.01, 

234 for which the photometric estimate of T^, is 459 K hotter than that of Mann et al., who find 

235 T^, ^ 3500K. The spectral type is MO. 5, which is marginally consistent with either of these 

236 values of T*. One explanation for the discrepancy is confusion with another star - there is 

237 a source about 4 arc-sec away in high-resolution UKIRT J-band imaging. The case of M 



238 dwarf stars is discussed in 



239 3.2. Habitable Planets 

240 Of the 2740 confirmed and candidate KOIs, the analysis of 1 star (KIC 7746948 hosting 

241 KOIs 326.01 and 326.02) failed, as it is missing an r magnitude and therefore cannot be 

242 analyzed by this procedure. Figure |5] is a Hertzsprung- Russell diagram of the other 2035 

243 KOI host stars. The stars hosting HZ candidates are indicated in black and only their 

244 associated errors in luminosity are plotted (for clarity). Seven HZ candidates (KOIs 375.01, 

245 422.01, 435.02, 490.02, 1096.01, 1206.01, and 1421.01) were excluded because their reported 

246 orbital periods have large uncertainties, being based on the duration of a single transit 

247 and the assumption of a circular orbit). The remaining HZ candidates have host stars 

24B with effective temperatures from 3500 K (early M-type) to 6100 K (late F-type). With the 

249 exception of KIC 6766634 (hosting KOI 1375.01), which is missing an r — z color, all HZ 

250 candidate host stars have all 6 photometric colors. Luminosities for some host stars have 

251 very high upper bounds because neither the photometric colors nor the priors can rule out 

252 the possibility that they are evolved with 95% confidence. But all are most likely to be 

253 dwarfs except for two of the hottest host stars (for KOIs 1168.01 and 1574.02), which have 
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254 probabilities of 77 and 53% of being evolved stars, i.e. logg < 4.2. All are designated as 

255 slighly metal poor, most likely the result of the metallicity prior ( §2.2p . 

256 The 52 candidates with pnz > 0.5 and most probable incident stellar irradiation in the 

257 HZ limits are listed in Tabled] and plotted in Figs . [6] an d [71 Thirty-four were previously 



258 identifi ed as possible HZ pla nets by 



Borucki et al. 



and/or 



Batalha et al. 



mm . iKaltenegger fc Sasselovl fl201lh 



2OI2I ). Most, but not all, of the others are candidates from the 



260 1-7-2012 release. The candidate around the brightest host star, KOI 87.01/A'e»/er-22b 

I 1 I 

261 was previousl y flagged bv iBorucki et al.1 (12011 ) and Kaltenegger fc Sasselov ( 2011 ) and 



262 confirmed by 



Borucki et al. 



(I2OI2I ). The photometric estimate of effective temperature 



263 (5735 K) is consistent with two spectroscopic estimates (5518 and 5642 K), the inferred 

264 luminosity is slightly higher 0.9Lq compared to 0.79Lq, and the inferred age of 8 Gyr 

265 i s consistent Avith sl ow rotation and low flux in the core of the Ca II H and K lines 

266 (IBorucki et al.ll2012l ). The inferred stellar mass is identical (O.93M0) that determined 

267 by astroseismology. The preferred planet radius is 2.32R^and is within the errors of the 

268 previously published value of 2.38 ± 0.13R®, although the 95% confidence interval for this 

269 star is large. 

270 KOI 250.04 is not (yet) a confirmed planet but is the outermost known member of 

271 the 4-pla net Kepler-26 syst em containing two components (b and c) confirmed by TTV 



272 analysis (ISteffen et al 



2OI2I ) and a fourth candidate (KOI 250.03) on the innermost orbit. 
273 The orbital period of KOI 250.04 {Pk = 46.83 d) is suspic ously close to one half of the 



period of a TTV signal seen near 90 d (ISteffen et al. 



2OI2I ). This analysis indicates that 



275 the host star of these planets has = 4072 K, i.e. is a late K dwarf. This is confirmed 

276 by two moderate-resolution visible-wavelength spectra which return 3996 K and 4067 K 



277 and a spectral 



type of K7.5 (Mann et al. in prep. 



T, = 3887K (IMuirhead et al. 



2OI2I) 



Steffen et al 



and an infrared spectrum which gives 



torn report T, = 4500 K based on 
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279 an SME analysis (jValenti &: Piskunovl Il996l ) of a Keck-HIRES spectrum. However, SME 

280 effective temperatures are unreliable for very cool stars such as this. KOI 250.04 has a 

281 radius of about 2.4R^, and it is of particular interest because further TTV analysis might 

282 constrain its mass. 

283 The distribution with planet radius among these 52 objects peaks in the super-Earth 

284 range (~ 2.5R®) and decreases with increasing radius, although there appears to be a 

285 second group of candidate planets with radii approximately that of Jupiter. Most of the 

286 smaller planets orbit lower- mass stars (Fig. [7j), a correlation at least partly due to detection 

287 bias. Presumably gas giants, these ob jects are potential hosts for habitable satellites 



288 (IKipping et al. 



2009 



Kaltenegger 



20101 ). The upper limits (95% confidece) or i the radii of 6 



289 candidates significantly exceeds the theoretical upper limit for cool Jupiters (IFortney et al. 



2010 



, Rp ~ 13. 5R®, ). Using the same approach as the calculation of pnz, I calculated 
the probability that a candidate planet's radius is larger than 1.2Rj, based on the transit 
depth and model stellar radii. For each of these 6 objects, there is a > 10% chance that the 
actual radius does exceed this limit. In 5 of these cases, this can be explained by the very 
large errors in the radius of the host star (an evolved star cannot be ruled out). The 95% 
confidence interval in radius of KOI 113.01 extends from 1.3Rj to a clearly non-planetary 
O.37i?0. The probability that this object has a radius larger t han 1.2-Rj is 9 9 %, su ggesting 
this is not a planet at all, but an very cool or ultracool dwarf. iBatalha et al.l ( 120121 ) list this 
KOI as having a "V-shaped" transit lightcurve indicative of an eclipsing binary. Depending 
on this object's luminosity, it is not inconceivable that any planets orbiting on stable orbits 
around it would be in the HZ of the combined system. All other KOIs in Table [H have 
p{Rp > l.2Rj) < 0.17. 

KO 3010.01 has a most probable radius of 1.35±^Q;o9Re; and an irradiance of ~0.81 
that received by Earth. Tabled] also includes 6 additional objects of interest with Rp < 211^) 
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304 planets and pnz > 0.1 (but < 0.5). Five of these orbit late K or early M-type dwarf 

305 stars, a figure that supports claims that these stars are the most promising locales to find 

306 Earth-size, and Earth-like planets. 



3.3. Not-so-habitable planets 



Kaltenegger fc Sasselovi ( 120111 ) list 27 planets with semimajor axes between the inner 
edge (as defined by the onset of a runaway greenhouse) and outer edge of the HZ. Of these. 



310 7 ( KOIs 113.01, 



the 



Batalha et al. 



6 5.01, 1008.01, 1026.01, 1134.02, 1168.01, 1232.01, were not retained in 



(120121 ) catalog. KOIs 1008.01 has a V-shaped transit shapes and 1232.01 
has a large radius indicative of an elipsing binary. KOI 1134.02 exhibits "active pixel offset" 
meaning that the target star is not the sourc e of the transit sign al. KOI 1026.01 might be 



an artifact of systematics in the Kepler data (iBatalha et al. 



20121). KOIs 11 3.01, 465.01 and 



Borucki et al. 



( 1201 ih catalog. Of the 



1168.01 were detected only with a single transit in the 
remaining 20, five (KOIs 139.01, 1099.01, 1423.01, 1439.01, and 1503.01) havepnz < 0.5 and 
so do not appear in this catalog, although KOI-1423.01 is omitted marginally only so (0.47). 
KOI 1439.01 is most strongly ruled out {phz = 0.06 ) because the revised i s 274 K hotter 



319 and R^: is 46% larger than the KIC values used by 



Kaltenegger fc Sasselovi (120111). The 



320 other 15 KOIs are retained in this catalog, along with 5 others from 



torn 



Kaltenegger &: Sasselov 



Batalha et al. 



(I2OI2I ) is problematic because 



322 A comparison with the HZ candidates of 

323 they use an equilibrium temperature criterion which is and dependent on the color /effective 
temperature of the host sta r. How ever, of the 24 candidate planets with 185K < Teq <300K 



325 in Table 8 of 



Batalha et al. 



mm . six (KOIs 119.02, 438.02, 986.02, 1938.01, 2020.01, and 

326 2290.01) have pnz < 0.05 and most probable < I > outside the HZ limits. In each case, 

327 this is because the new estimates for T^, are >200 K hotter than the previously published 
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32B values, and becuse the most probably estimate of radius is significantly larger. 



329 4. Discussion 

330 Assumptions and biases: My priors do not account for the probability that a planet 

331 will transit its host star and be detected by Kepler, and thus be included in the KOI 

332 catalog. S uch selection effects can be important in catalogs of transiting planets and their 



333 host stars ( IGaidos &: Manrul2013[ ). The geometric transit probability R^/a, where a is the 

334 semimajor axis, is proportional to r^/^ and could be included readily enough: this factor 

335 will favor stellar models with larger radii. However, the probability of transit detection 

336 is primarily related to transit depth 6 ~ {Rp/R^fY and for a given 5, a prior on stellar 

337 radius is ultimately a prior on planet radius. Some of these KOIs are nearly Earth-size, 

338 where the completeness of the Kepler survey is still being refined. Other KOIs are at or 

339 near theoretical limits of giant planet radii and any prior on stellar radii would have to 

340 include scenarios for astrophysical false positives. There are additional, but perhaps minor 

341 complexities: the probability of a transit occurring and being detected will also depend 

342 on e, 0, 6, as well as D, the transit duration. For these reasons, I do not include transit 

343 detection as a prior. 

344 The derivation of the extinction (Eqn. [2]) presumes a linear relationship between 

345 extinction in different bandpasses, i.e. that all can be linearly related to reddening Eb-v- 

346 This is not strictly correct, but is a fair approximation in the limit of small reddening. 

347 The median derived Eb-v for these stars is only 0.069, corresponding to 0.2 magnitudes 
34B of extinction, and the 95 percentile value is 0.18. The distribution of reddening values is 

349 shown in Fig. |H1 The scale height of dust is smaller than that of stars, and as a result the 

350 assumption of identical gas and dust distributions (Appendix A) slightly underestimates 

351 the amount of reddening, as well as the temperature and luminosities of Kepler stars. 
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352 The total number of candidate HZ planets is not sensitive to the precise irradiation 

353 limits. Because of detection bias towards short-period orbits, there are essentially no 

354 detected planets beyond the HZ (Fig. [H]), except perhaps for those detected by only a single 

355 transit. For a planet with 100% cloud cover, the runaway greenhouse irradiation limit is 

356 23% higher than the 50% case, but this admits only one additional candidate to the catalog. 

357 The trouble with M dwarfs: The fundamental parameters of M dwarf stars have been 

358 a notorious challenge for models because of the difficulty in reproducing the observed 

359 mass-radius relation and their complex spectra. The DSED models employed here 



accurately predict s the radii o: 
system KOI-126 ( IFeiden et al. 



362 ( iHauschildt et al. 



the t wo M dwarfs in the triply eclipsing hierarchical triple 



201l|). DSED uses the PHOENIX model atmospheres 



19991 ) for both the stellar surface boundary conditions and to generate 
synthetic magnitudes. The analysis of Mann et al., in prep., compares the BT-SETTL flavor 
of PHOENIX models to medium-resolution spectra, hence the consistency between effective 
temperatures estimates is not entirely surprising. However, not all flavors of PHOENIX are 
identical, and moreov er, Mann et al. calibrated temperatures against the stars analyzed in 



Boyajian et al. 



( I2OI2I ). which have both measured parallaxes and (interferometric) radii. 

36B For early M dwarfs, a difference of only ~100 K in T^, is nevertheless a difference in of 

369 about 30%. 

370 Another obstacle is that accurate modeling of the lightcurves of planets transiting 

371 M dwarfs must correctly account for significant limb darkening in the Kepler pass-band. 

372 Erroneous transit durations, acting through the prior described in §2.3^ can bias the analysis 

373 towards models with erroneous radii: a 10% error in R^, again leads to a ~30% error in L^,. 

374 This is sufficient to move a planet completely outside the HZ, or at least decrease the pnz 

375 of a marginal HZ planet to <50%. A re-analysis of the Kepler transit lightcurves with an 

376 improved limb-darkening model and re-derivation of the parameters of M dwarf KOI hosts 
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377 is forthcoming (Mann et al., in prep.) 



Stellar parameters based on analysis of photometry are no substitute for values 



379 based on high-reso 



Pinsonneault et al. 



ution spectra, as long as the latter are carefully calibrated (e.g., 



2OI2I ). However, the median magnitude of the host stars of these planet 



381 is Kp 15.1, and high-resolution spectroscopy would be very observationally expensive. 

382 The object most amenable to foUowup is, not coincidentally, Kepler-22h {Kp = 11.7). The 

383 next brightest host star is that of KOI 371.01 {Kp = 12.4) and the rest are much fainter still 

384 and would require significant time on very large telescopes. The analysis presented here 

385 does, however, generate a robustly-defined catalog of targets to prioritize such work. 

386 The Gaia (originally Global Astrometric Interferometer for Astrophysics) mission, 

387 scheduled for launch in October 2013, will obtain parallaxes with a sky- averaged, end-of- 

388 mission precision of 25 yuas and 40 /ias for 15th and 16th magnitude stars, respectively. 



with somew! 



390 ( Ide Bruijne 



lat s uperior performance at the ecliptic latitude (~66 deg.) of the Kepler field 



2OI2I ). To assess the potential of Gaia to refine the habitable zones of Kepler 
stars and the sizes of the planest that inhabit them, I re-calculated Pj (Eqn. [T]) for all 
models using a prior for distance modulus based on Gam's expected precision: 

(/i - /io)^ 



Pf, = exp 



11 



where /iq is the most probable distance modulus from the original analysis, and 
(T^ ~ 8.7 X 10'^''/^"^ is the uncertainty in distance modulus from a 40 /xas precision in 
parallax. The 95% confidence intervals in radius and stellar irradiance of the 52 HZ 
candidates were re-calculated and are shown in Fig. [91 The most probable values are 
not changed, but the fractional errors in radius and irradiance are reduced by a factor 
of ~5, from a median of 10% and 24%, respectively, to 1.7% and 5%, (equating 95% 
confid ence intervals to 4cr). The largest planets tend to orbit the hottest and most distant 



stars (iGaidos fc Mann 



20131) and their parameters have the largest errors. Typically, a 
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401 few hundred models contribute to the calculation, i.e. retain appreciable P values, but in 

402 a few cases the number is a few dozen and finite model grid size may contribute to the 

403 error budget. Values of pnz for 50 of the 52 planets are > 0.97. Spectroscopic values of 

404 T^, accurate to 100 K would offer only modest further improvement (1.5% and 4.5% errors, 

405 respectively). Because these precisions reach or exceed levels of confidence in the predictions 

406 of the stellar models themselves as well as the absolute calibration of the photometry, 

407 refinement and verification of these may prove a more cost-efficient avenue for improvement. 



For example, Ugri and p hotometry o 



the Isaac Newton Telescope ([Greiss et al 



at WIYN flEverett et al. 



2012f ). 



much of the Kepler field has been obtained at 
2OI2I ) and UBV photometry has been obtained 



411 With such precision (or better!), it should be possible to locate planets within different 

412 realms of the HZ, i.e. near the inner edge, where low CO2 atmospheres, and possibly high 



413 cloud fraction if t 

414 where high CO2 ( 



lere is a temperature -cloud feedback, should prevail: or the outer edge. 



von Paris et al. 



20131 ). and possible water cloud-free atmospheres are 



415 more likely Candidate HZ planets in multi-planet systems might be confirmed or even have 

416 masses determined by TTV. Although such advances may be challenges for planets around 



or surveys of nearb y. 



417 faint Kepler stars, this analysis offers a preview of the potential 

418 more observationally accessible stars, e.g. by the proposed TESS (iDeming et al. 

419 CHEOPS missions. 



200% and 
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492 A. Derivation of a uniform probability distribution for extinction 

493 If the probability distribution of stars with distance x along the line of sight is f{x) 

494 and the density of dust is g{x), then the total column density of dust along the line of sight 

495 to a particular star is 



496 where Aq is a constant factor. The probability of extinction to any randomly selected star 

497 falling between A and A + dA is 



498 However, from Eqn. lAll dx/dA is simply g{x)~^ and if f{x) and g{x) are identically 

499 distributed with x, then p{A) is a constant, i.e. uniformly distributed over the range of 

500 allowed values. 
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Fig. 1. — Priors of age, metallicity, and distance modulus generated from 428,792 
TRILEGAL-simulated stars with Kp < 16 and logg' > 4 in the Kepler field. The distri- 
butions of all synthetic stars is presented here, but only the ~10% of stars with Galactic 
latitudes within 0.5 deg of a Kepler star of interest are used to generate actual priors. 
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1.2 




Fig. 2. — Priors on the planet transit duration as a function of the parameter A = 
D/{t^/^P]1^), where D is the transit duration, r is the stellar free-fall time, and Pk is the 
Keplerian orbital period, for circular orbits (solid line) and orbits with Rayleigh-distributed 
eccentricities with means of 0.1, 0.2, and 0.3. 
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Fig. 3. — Comparison between effective te mperatures based on ph otome try and spectroscopi c 



values. Circles are solar-type stars from 



Buchhave et al. 



(120121 ) and 



Bruntt et al. 



(120121). 



Triangles are M dwarfs from Mann et al., in prep.. Black points are stars where all six 
photometric colors are available; grey points represent stars where at least one color is 
unavailable. The solid line i s equality between t he est imates and the dashed line represents 



the ~215 K offset found by 



Pinsonneault et al. 



torn . 
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Fig. 4. — Comparison between photometric an d spec trosc opic estimate s of su rface gravities 



where the latter are taken from 



Buchhave et al. 



(120121) and 



Bruntt et al 



( 120121 ). Black points 



are stars where all six photometric colors are avialable; grey points represent stars where at 
least one color is unavailable. The solid line is equality between the two estimates. 
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Fig. 5. — Hertzsprung- Russell diagram for host stars of 2739 candidate and confirmed Kepler 
planets. Black points are the 52 candidate HZ planets in Table 1. The error bars represent 
95% confidence intervals on luminosity and are shown only for the candidate HZ planets. 
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Fig. 6. — Radius and stellar irradiance of candidate and confirmed Kepler planets, and the 
Earth. Candidate planets in habitable zones are highlighted as black, all other KOIs are 
grey, and a the vast majority of KOIs experience higher irradiances and fall outside the 
left-hand boundary of the plot. The error bars correspond to 95% confidence intervals. The 
solid, dotted, and dashed hues are the boundaries of the HZ for a 50% cloud-covered Earth 



like planet aro und a solar 



(7800 K) star flSelsis et al 



type 



20071 ). 



star (5780 K), an early M dwarf (3800 K) and late A-type 
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Fig. 7. — Luminosity of the host star vs. orbital period of candidate HZ planets detected 
by Kepler^ plus the Earth. The points are scaled to planet radius and the the darker the 
point, the more likely it is in the HZ. T he two lines delim it the boundaries of the HZ for 



Earth-like planets with 50% cloud cover (jSelsis et al 



20071 ). To plot these boundaries with 



these axes, it was necessary to assume simple but standard power-law relations between the 
luminosities, masses, and effective temperatures of main-sequence stars. 




Fig. 8. — B — V interstellar reddening along the line of sight to Kepler KOI host stars 
calculated using Eqn. |2]for the best-fit stellar models. The median value is 0.069. 
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Fig. 9. — Reduction in the uncertainties in planet radius and stellar irradiation expected 
from inclusion of Gaia parallax measurements with 40 /xas errors. Compare to Fig. [61 



Table 1. Candidate Planets in the Habitable Zones of Kepler Stars 



Planet Parameters Stellar Parameters 



KOI 


KIC 


Phz 


Period 


Irrad. 


Planet Radius (R®) 




log g 


[Fe/H] 






Age 


Comment'' 








(days) 


(/ ) 


All'" 


lo\v(>r 


U]>|)('1' 


(K) 






[L . ) 


(.M-) 


(Gyr) 




87.01 


10593626 


0.64 


289.86 


1.29 


2.32 


2.00 


5.91 


5735 


4.43 


-0.2 


0.90 


0.93 


7.8 


Boll,KSll,Bul2,Xep«er 


113.01 


2306756 


0.76 


386.59 


0.62 


14.35 


14.33 


40.54 


5575 


4.56 


-0.2 


0.62 


0.89 


1.7 


KSll,V-shaped transit 


250.04 


9757613 


0.98 


46.83 


1.72 


2.40 


1.97 


2.54 


4072 


4.71 


-0.3 


0.07 


0.55 


3.7 


Mal3,/<:ep/er--26 


351.01 


11442793 


0.89 


331.64 


1.54 


11.10 


9.63 


13.91 


6012 


4.34 


-0.3 


1.44 


1.10 


6.1 


Boll.KSll 


371.01 


5652983 


0.70 


498.39 


0.72 


4.21 


3.67 


19.24 


6001 


4.48 


-0.1 


1.09 


1.00 


1.5 




401.02 


3217264 


0.54 


160.02 


1.96 


3.76 


3.30 


4.12 


5454 


4.50 


-0.2 


0.59 


0.87 


7.5 


Boll,KSll 


433.02 


10937029 


0.99 


328.24 


0.87 


11.11 


9.58 


12.52 


5650 


4.52 


-0.2 


0.71 


0.92 


3.7 


Boll.KSll 


518.03 


8017703 


1.00 


247.35 


0.58 


2.82 


2.52 


3.04 


5030 


4.63 


-0.3 


0.28 


0.74 


1.4 




622.01 


12417486 


0.77 


155.04 


1.65 


6.01 


5.22 


23.00 


5256 


4.53 


-0.1 


0.46 


0.83 


7.3 


Boll.KSll 


682.01 


7619236 


0.93 


562.14 


0.74 


8.56 


6.97 


17.11 


6114 


4.43 


-0.1 


1.40 


1.10 


1.4 


KSll,Bul2 


701.03 


9002278 


1.00 


122.39 


1.52 


2.04 


1.80 


2.16 


5026 


4.63 


-0.3 


0.29 


0.72 


1.2 


Boll,KSll,Bul2 


812.03 


4139816 


0.65 


46.18 


1.87 


2.34 


1.95 


2.38 


4084 


4.69 


-0.3 


0.08 


0.57 


4.0 


Boll,KSll,Mul2,Mal3 


854.01 


6435936 


1.00 


56.06 


0.70 


2.10 


1.95 


2.42 


3658 


4.79 


-0.1 


0.04 


0.48 


1.4 


Boll,KSll,Mul2,Mal3 


868.01 


6867155 


0.78 


235.97 


0.37 


10.77 


8.16 


11.13 


4474 


4.64 


0.2 


0.17 


0.74 


1.4 


KSll,Mul2,Mal3 


881.02 


7373451 


0.99 


226.89 


0.88 


4.39 


3.95 


4.89 


5225 


4.57 


-0.3 


0.40 


0.78 


6.5 


KSll 


902.01 


8018547 


1.00 


83.91 


1.72 


6.51 


5.97 


7.06 


4623 


4.63 


-0.2 


0.19 


0.68 


6.1 


Boll,KSll,Mul2 


1168.01 


10460629 


0.61 


432.13 


1.37 


4.18 


3.47 


6.30 


5855 


4.15 


-0.1 


2.09 


1.35 


8.5 


KSll 


1209.01 


3534076 


0.78 


272.07 


1.16 


7.06 


5.98 


18.26 


5668 


4.50 


-0.3 


0.72 


0.89 


5.5 


Bal2 


1298.02 


10604335 


1.00 


92.74 


0.93 


2.14 


1.76 


2.38 


4262 


4.67 


-0.3 


0.11 


0.60 


6.3 


Mal3 


1356.01 


7363829 


0.92 


384.03 


0.98 


8.81 


7.61 


14.34 


5809 


4.41 


-0.2 


1.05 


1.00 


6.2 




1361.01 


6960913 


0.66 


59.88 


1.23 


2.14 


2.03 


2.89 


4140 


4.73 


-0.5 


0.07 


0.53 


2.8 


Boll,KSll,Mul2,Mal3 


1375.01 


6766634 


0.67 


321.21 


1.29 


5.13 


4.45 


10.83 


6001 


4.48 


-0.1 


1.09 


1.00 


1.5 


Boll.KSll 



Table 1 — Continued 



Planet Parameters Stellar Parameters 



KOI 


KIC 


Phz 


Period 


Irrad. 


Planet Radius (R®) 




log g 


[Fe/H] 




M, 


Age 


Comment'' 








(days) 


{%) 


MP" 


lower 


upper 


(K) 








{Me) 


(Gyr) 




1429.01 


11030711 


0.77 


205.92 


1.79 


4.98 


4.26 


7.30 


5623 


4.45 


-0.2 


0.80 


0.94 


8.0 


Boll,KSll 


1430.03 


11176127 


0.88 


77.48 


1.58 


2.67 


2.07 


2.87 


4532 


4.66 


-0.3 


0.16 


0.69 


1.2 


Bal2 


1431.01 


11075279 


0.88 


345.16 


0.72 


6.75 


6.20 


15.41 


5477 


4.50 


-0.2 


0.61 


0.87 


7.5 


Bal2 


1466.01 


9512981 


0.96 


281.56 


0.51 


12.04 


9.68 


13.33 


4892 


4.55 


0.0 


0.30 


0.76 


11.8 


Bal2 


1477.01 


7811397 


0.97 


339.08 


0.52 


9.70 


8.47 


11.68 


5236 


4.56 


-0.2 


0.42 


0.82 


6.0 


KSll 


1574.02 


10028792 


0.99 


573.96 


1.12 


5.22 


3.99 


6.23 


6131 


4.19 


-0.2 


2.45 


1.32 


6.0 




1582.01 


4918309 


0.67 


186.40 


1.92 


6.19 


5.14 


14.62 


5668 


4.50 


-0.3 


0.72 


0.89 


5.5 


Boll,KSll 


1596.02 


10027323 


0.93 


105.35 


1.39 


2.78 


2.36 


2.94 


4666 


4.62 


0.2 


0.22 


0.77 


1.8 


Boll,KSll,Mul2 


1686.01 


6149553 


1.00 


56.87 


0.56 


1.16 


0.86 


1.44 


3598 


4.84 


-0.2 


0.03 


0.44 


1.8 


Bal2,Mal3 


1739.01 


7199906 


0.62 


220.66 


1.47 


1.96 


1.74 


4.13 


5650 


4.52 


-0.2 


0.71 


0.92 


3.7 


Bal2 


1871.01 


9758089 


1.00 


92.72 


1.48 


2.46 


2.31 


2.65 


4649 


4.64 


-0.2 


0.19 


0.71 


1.4 


Bal2 


1876.01 


11622600 


0.99 


82.53 


1.66 


o.lo 


O c^7 
Z.O t 


o.ol 


4622 


4.64 


-0.2 


0.18 


0.72 


2.1 


Bal2 


1902.01 


5809954 


0.91 


137.87 


0.36 


2.63 


1.93 


2.63 


3886 


4.71 


0.2 


0.07 


0.58 


1.0 


Bal2,Mal3 


2102.01 


7008211 


0.83 


187.75 


0.92 


3.01 


2.94 


12.37 


5165 


4.62 


-0.5 


0.29 


0.68 


7.7 


Bal2 


2418.01 


10027247 


1.00 


86.83 


0.63 


1.62 


1.31 


1.85 


3959 


4.74 


-0.4 


0.06 


0.53 


1.2 


Bal2,Mal3 


2469.01 


6149910 


0.55 


131.19 


0.99 


2.16 


2.17 


8.92 


4704 


4.64 


-0.2 


0.21 


0.73 


1.6 


Bal2 


2626.01 


11768142 


1.00 


38.10 


0.72 


1.14 


0.97 


1.61 


3514 


4.90 


-0.3 


0.02 


0.38 


1.5 


Bal2,Mal3 


2681.01 


6878240 


0.62 


135.49 


1.33 


5.38 


5.01 


8.00 


5026 


4.63 


-0.3 


0.29 


0.72 


1.2 




2686.01 


7826659 


0.84 


211.04 


0.50 


3.24 


3.12 


21.08 


4647 


4.64 


-0.2 


0.19 


0.70 


1.5 




2689.01 


10265602 


0.56 


165.35 


1.95 


6.04 


5.01 


14.01 


5544 


4.53 


-0.1 


0.64 


0.92 


3.5 




2691.01 


4552729 


0.91 


97.46 


1.82 


3.63 


3.36 


12.11 


4842 


4.62 


-0.2 


0.26 


0.75 


1.6 




2703.01 


5871985 


0.99 


213.26 


0.48 


3.49 


3.18 


3.69 


4549 


4.63 


0.2 


0.19 


0.71 


1.1 





Table 1 — Continued 



Planet Parameters Stellar Parameters 

KOI KIC PHZ Period Irrad. Planet Radius (R®) T* log g [Fe/H] L* M* Age Comment'' 

(days) (/gj) MP" lower upper (K) (Lq) (Mq) (Gyr) 



2757.01 


6432345 


0.75 


234.63 


1.41 


2.92 


2.57 


6.21 


5668 


4.50 


-0.3 


0.72 


0.89 


5.5 




2762.01 


8210018 


1.00 


132.99 


0.92 


2.50 


2.33 


2.70 


4647 


4.64 


-0.2 


0.19 


0.71 


1.6 




2770.01 


10917043 


1.00 


205.38 


0.37 


2.33 


2.16 


2.70 


4368 


4.68 


-0.3 


0.12 


0.62 


1.2 


Bal2 


2834.01 


5609593 


0.96 


136.21 


1.04 


2.67 


2.35 


2.94 


4764 


4.63 


-0.1 


0.23 


0.74 


1.5 




2933.01 


12416987 


1.00 


119.08 


1.05 


3.84 


2.72 


3.99 


4549 


4.63 


0.2 


0.19 


0.71 


1.1 




2992.01 


8509442 


1.00 


82.66 


0.57 


2.11 


1.94 


2.81 


3923 


4.77 


-0.4 


0.05 


0.50 


2.0 




3010.01 


3642335 


1.00 


60.87 


0.81 


1.35 


1.26 


1.79 


3918 


4.79 


-0.5 


0.05 


0.49 


2.0 




3086.01 


10749059 


0.84 


174.74 


1.81 


3.27 


2.76 


5.70 


5537 


4.53 


-0.1 


0.64 


0.91 


3.3 














Other Planets with R 


p < 2R0 and pnz > 0.1 












775.03 


11754553 


0.13 


36.45 


2.50 


1.93 


1.73 


2.16 


4132 


4.71 


-0.4 


0.08 


0.56 


3.6 


Bul2, 


1422.02 


11497958 


0.39 


19.85 


2.39 


1.85 


1.34 


2.01 


3638 


4.83 


-0.3 


0.03 


0.45 


1.4 




2124.01 


11462341 


0.36 


42.34 


2.42 


1.19 


1.05 


1.41 


4229 


4.70 


-0.4 


0.10 


0.59 


1.5 


Bal2,Mal3 


2650.01 


8890150 


0.48 


34.99 


2.77 


1.61 


1.12 


1.63 


4039 


4.67 


0.1 


0.09 


0.61 


5.5 


Bal2 


2931.01 


8611257 


0.20 


99.25 


2.31 


1.89 


1.68 


7.60 


5126 


4.61 


-0.2 


0.35 


0.79 


1.5 




3034.01 


2973386 


0.46 


31.02 


2.11 


1.61 


1.40 


2.01 


3923 


4.77 


-0.4 


0.05 


0.50 


2.0 





^MP = most probable radius 
''Reported as HZ candidate in: Bll 



reported in Bul2 



Buchhave et al 



Borucki et al. 



I201J), Mul2 



I 2011h. KS12 



Muirhead et al 



Kaltenegger fc Sasselovl J2011h . Bal2 



Batalha et al. 



12012); Spectroscopy 



(201J), Mal3 = Mann et al., in prep. 



